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ABSTRACT. It is shown that, if the vorticity magnitude associated with a (presumed singular) three-dimensional incom- 
pressible Navier-Stokes flow blows-up in a manner exhibiting certain time dependent local structure, then time independent 
ff) i estimates on the L 1 norm of \u\ log y/1 + \uj\ 2 follow. The implication is that the volume of the region of high vorticity 

decays at a rate of greater order than a rate connected to the critical scaling of one-dimensional local sparseness and, conse- 
C ~_~) ' quently, the solution becomes sub-critical. 

(N ' 



I. Introduction 

The purpose of this article is to present blow-up scenarios under which Leray solutions to the 3D Navier-Stokes 
, equations (3D NSE) behave sub-critically with respect to the critical scaling of one-dimensional local sparseness. In order 
(— | ' to contextualize our discussion we begin with several general remarks. Speaking informally, there is a scaling-gap 
between known a priori finite quantities - e.g. sup 0<t<T | |u| \\ and sup 0<t<T | |Vu||i - and quantities on which regu- 
larity can be conditioned - e.g. sup 0<t<T | |u| I3 and sup 0<t<T 1 1 Vw| ||. Under the natural scaling for 3D NSE, 

, . 1 fx t 

>' 

t — we see that the a priori controlled quantities exhibit sub-critical scaling - our examples scale as A - while those 
sufficient for regularity scale critically as A . This mismatch is referred to as the "scaling-gap" and indicates the 
■ Navier-Stokes problem is a super-critical problem. 

£C) ' A theme apparent in many regularity results is the inclusion of premises that explicitly bridge the scaling gap. For 
example, a significant result in regularity theory was non-existence of (backward-in-time) self -similar blow-up (cf . 
l|T8ll23l[T4l ). Because self -similar solutions are scaling invariant in virtue of their construction, their study is effec- 
tively a restriction to a class of critical solutions. Indeed, self -similar solutions satisfy a scaling invariant point- wise 
bound, esssup s ((|a; — xq\ + y/T — t)\u(x, t)\) < 00. Regularity of solutions satisfying this estimate has been af- 

k>( firmed if the solution is additionally assumed to be axisymmetric (cf. Bll3l l20l[T6ll ) but the general case remains an 
5_J open problem. Regarding the axisymmetric case, the theme originally referenced is again apparent: first, the flow 
is assumed to exhibit some feature which is critical and, second, additional restrictions are identified from which 
regularity follows. Similar examples (in particular those from the e-regularity theory) require additional conditions 
on the smallness of the presumed finite scale-invariant quantity. 



In two recent complementary publications, [12 8 ], a new dynamic approach to bridging the scaling gap is elucidated. 
The physical motivation is the persistence (in the average sense) of the axial lengths of vortex filaments. This picture 
is supported in experimental and numerical studies and by a mathematical result which is the main result of El. 
By considering this length persistence and the decay rate for the volume of the region of intense vorticity, a connection 
is found between the scaling of the latter quantity and the critical scaling of one-dimensional local sparseness. As this 
argument provides the context for our own results we make the matter more precise. 

The (rigorous) regularity criteria is presented in [12 J. In that paper, regularity of mild solutions to 3D NSE with 
L°° (M 3 ) initial data is conditioned on a geometric measure-type criteria involving one-dimensional linear sparse- 
ness of the super-level sets of the vorticity. The implication is that the critical scaling for the local one-dimensional 
sparseness of the region of intense vorticity is of order Co ||w||oo ■ The proof of this geometric measure-type regular- 
ity criteria is based on an interplay among the diffusion, the basic symmetries present in the 3D NSE and geometric 
properties of the harmonic measure (resulting in strong anisotropic diffusion); this in turn draws on the ideas ex- 
posed in ||TT|. On the other hand, the a priori estimates on the L 1 norm of the vorticity found in JHIlZI imply the 
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volume of the region where vorticity magnitudes are high decays according to, 

Vol(A,(>( t , IU ))<K 

where A t (y) = {x : \u>(x, t)\ > y}. Interestingly, if the region of intense vorticity corresponds precisely to the space 
occupied by filamentary vortex structures and the length of these filaments is non-decreasing - that is, pinned to 
the characteristic length scale of a turbulent region - then the anti-axial diameters of these filaments exhibit a rate of 

— 1/2 

decay of order at least C3 1 \uj(t) \ \oo , which matches the critical scaling for local one-dimensional sparseness. It is in this 
sense that the problem is rendered critical. That the vortex filaments have persistent lengths (in the average) which 
are comparable to the scale of the turbulent region is not yet rigorously established but is supported by numerical 
evidence as well as a mathematical evidence which is the prime consideration of JSJ. 

Even assuming the soundness of the above argument it is not immediately evident that regularity follows (it would 
if one could reconcile the possible difference between co and C3). The present paper partially overcomes this by 
illustrating two blow-up scenarios in which the region of intense vorticity decays at a faster rate than the critical rate 
discussed in [12. 8], thereby rendering the constants irrelevant. In particular, we are interested in concluding, 

voi(A t (l|Mi)iu)) < c ; 

where, in our first scenario, <£>(.t) = log(l + x), and, in the second, = loglog(e + x). These results are given 

in Section [3] We note for clarity that ours are not stated as regularity criteria but instead scenarios under which 3D 
NSE becomes sub-critical with regard to the scaling described in |Tl2l l8l. 

The above decay rates will be obtained by imposing certain structural requirements on the blow-up rate(s) exhibited 
by vorticity components. Let lo^ denote the positive and negative vorticity components truncated away from zero. 
We will define amenable blow-up rates of orders and 1 in detail in Definition [6] and here only illustrate a class of 
functions exhibiting such blow-up profiles and describe how the blow-up rates are connected to . A function g 
(to be identified with one of the ^fs) exhibits a local algebraic blow-up (around xq, at time T) if there exists a constant 
C > 1 so that, for (x, t) in a cylinder, Q = B(x , r) x (0, T], we have, 

1 / 1 \ Q ' ±(t) / 1 

C{\pf{ x ,t)\+T±(t)) ~ l9{X ' t)l ~ C {\ P f(x,t)\+T?(t) 

where, at each time t, pf{-,t) is some polynomial of dimension less than a fixed natural number , orf is positive 
valued and bounded away from both and 00, and tj 1 is a positive (up to T) scalar function of time which vanishes 
at the singular time. There is a considerable amount of freedom present in the above construction as the polynomial 
is allowed to vary (wildly) in the time dimension. We will also consider a scenario where some asymmetry is al- 
lowed between the bounds assumed on vorticity components. In particular, for appropriate blow-up rates (x, t), 
we will require, 

Df{x,t) < \wf(x,t)\ < CDf{x,tft®, 

where fif : (0, T) — > [1, Bj < 00]. 

A key role is played in our analysis by the cancelations evident in the vortex stretching term in the context of the real 
Hardy space H 1 exploited via the Div-Curl lemma |5| and the H 1 — BMO duality flOl 9]. In the standard way this 
gives uniform-in-time control of the vortex stretching term. The structural blow-up assumptions are provided to 
ensure uniform-in-time control of the BMO norm of a multiplier (multiplied against the vortex stretching term). In 
the algebraic case, this is enabled by the logarithm's depletive effect on the unboundedness of the mean oscillations 
of polynomial functions. The effectiveness is witnessed by the remarkable fact (cf. ||T9l |2H ) that there exists a 
constant C(d, n) so that, for any polynomial on K™ of degree less than or equal to d, 

\\ lo S\P\\\ B MO^ C ( d ' n )- 

In particular, the constant is independent of the coefficients. This allows us to introduce time-dependent algebraic com- 
parability conditions on the spatial profiles prior to a possibly singular time and do so in a manner that preserve 
time-independent estimates on the BM O norms of the logarithms of these profiles. 
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We proceed in Section [2] to review needed results from harmonic analysis and then define the classes of blow-up 
scenarios which will be amenable to our PDE argument. We also include results which connect the structural blow- 
up assumptions to an energy inequality-type argument given in Section [3] The statements and proofs of the main 
results are contained in Section[3] 



2. Preliminaries and Amenable Blow-Up Rates 



Here we review needed results from harmonic analysis and present a lemma which will connect these ideas to the 
PDE context of Section|3] The maximal function of a distribution / is defined for all x G R" as, 

M h f(x) =sup\f*ht(x)\, 
t>o 

where h is a fixed test function supported on the unit ball so that J h dx = 1. 
Definition 1. The distribution / is in the Hardy space W. 1 if 1 1/| |^i := | \Mhf\ | i < oo. 

In J5)/ Coifman, Lions, Meyer, and Semmes reformulated and refined some key features of the 'sequential' theory of 
compensated compactness within the framework of Hardy spaces, the key idea being that certain nonlinear quantities 
exhibit cancelations yielding the improved regularity. One such result is the Div-Curl lemma. 

Lemma 2. (Coifmann, Lions, Meyer, Semmes - (51) Suppose E,B e L 2 (R 3 ) 3 with V • E = and V x B = (in the sense 
of distributions). Then, E ■ B etf, and there exists a universal constant C such that 

\\E-B\\ H i < C\\E\\ L 2\\B\\ L 2. 

Because weak solutions to the Navier-Stokes equations are divergence free and because the curl of a gradient is 
always equal to the trivial distribution, it follows that the advective term, (u ■ V)it, in the velocity-pressure formu- 
lation exhibits div-curl structure, as does the vortex stretching term, (w ■ V)u, in the vorticity-velocity formulation. 
Consequently, using Hardy spaces, refined regularity results can be established for weak solutions of 3D NSE (see 
Chapter 3.2 of fT7\ for a collection of such results). 

The well known result of Fefferman (cf. Ifl0l l9l, as well as Stein's monograph [22j) establishes that the dual space 
of W 1 is precisely the space of functions of bounded mean oscillation, BMO. By f B we denote the average of the 
locally integrable function / over the ball B of volume \B\; i.e., f B = \B\ _1 J B f dx. The mean oscillation of / over 
B is l-Bl" 1 / |/ — f B \ dx and is the quantity typically used to characterize BMO. 

Definition 3. The locally integrable function / is in BMO if 1 1/| \bmo < °o where 

WIWbmo ■= sup -— / \f(y) - f B ( x ,r)\ dy. 

igS 3 ;0<r \- a ( x i r )\ JB(x,r) 

Several results about the space BMO will expedite future work. An elementary fact about this space is that, if 
g G BMO, then the following implication holds, 

(1) g(x) + m < f(x) < g(x) + A/ for a. e. x G M 3 ==> \\f\\ B MO < I Iff 1 1 bmo + M — m. 

The above is easy to see by directly comparing |/ — fs\ to \g + M — (g + m)s| and concluding that, for all balls B, 
we have, 

If, | 1 f , . 

\f - jB\dx< — I \g- g B \dx + M - m. 



\B\J B ^ J J ^ ~\B\ JB 

We include a technical lemma regarding certain point-wise multipliers on BMO which appears in a paper by 
Iwaniec and Verde (cf. Lemma 2.1 of Ifl5l ; see also Lemma 5.10 of [2]). This result will allow us to use the standard 
iV^-BMO duality in an essentially local context instead of working in the local versions of these spaces. It is worth 
mentioning that our results can be achieved in the local spaces using a local non-homogeneous version of the Div- 
Curl lemma (cf . J5l) paired with various local h 1 — bmo dualities (cf . [2]) thereby obtaining estimates on the localized 
vortex stretching term (this approach is similar to the estimates in 1|T3"10 . 

Lemma 4. (Iwaniec and Verde - |15p Suppose h G BMO(R n ) and </> 6 C^(B(x ,R)). Then cf>h e BMO(W l ) and we 
have, 



bmo <C{n,R)\\V4>\\c 



BMO 



mS B hdx 
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The space BMO enjoys an intimate connection with the logarithm. This is illustrated in the following lemma due 
to Stein. 

Lemma 5. (Ricci and Stein - fl9\ ) Let P be any polynomial in W 1 of degree less than or equal to d. Then there exists a 
constant C = C(d, n) so that log \P\ £ BMO and 

||log|P||sMO < C{d,n). 

The fact that the constant appearing above is independent of the coefficients evidences the remarkable extent to which 
the logarithm depletes mean oscillations and motivates the assumptions on which the conclusions of Section [3] are 
conditioned. This will be done by assuming that solutions blow-up in a fashion possessing certain structure. The 
remainder of this section is dedicated to describing amenable blow-up structures and the statement and proof of a 
lemma which will connect this structure to the PDE context of [3] We denote by log" 1 (2) an iterated composition of 
the logarithm with itself m — 1 times, that is, log™ (a;) = log o • • ■ o log(V e m + x' 2 ) where m applications the logarithm 
are carried out and e m is a power of e defined so that log m (0) = 0. We will use the following notation when defining 
functions in (x, t) possessing singularities at time T in the spatial set S C R 3 , 

Q T {x ,R, S) = B(x ,R) x (0, T) U ((B(x ,R) \S)x {T} and Cl T (S) = R 3 x (0, T) U (M 3 \ S) x {T}, 

where R > and xq £ M 3 . In contexts where xo and R are fixed we will make the abbreviation Qt(S) = 
Q T (x ,R, S). 

Definition 6. Fix x e K 3 , R, T > 0, m £ No and a set of measure zero, S, which is contained in a compact subset of 

B(xo,R). 

a. An amenable blow-up rate of order m on Qt(S) is a function, D(x,t) : £It{S) — > [0,oo), which additionally 
satisfies, 

i. there exists M > so that sup 0<(<T || log m+1 D(x, t)\\ B MO < M)/ 

ii. there exist M U M 2 > so that 1/Mi < D(x, t) on Q T (S) and D(x, t) < M 2 on B(x , R) c x (0, T\. 

b. A function f(x, t) : Qt(S) — > R exhibits an amenable blow-up rate of order on Qt(S) if there exists an amenable 
blow-up rate of order 0, D(x, t), on Qt(S), so that, for some C* > 1, 

D(x, t) < \f(x, t)\ < C*D(x, t) for all (x, t) G Q T (S). 

c. A function fix, t) : Qt(S) —> R exhibits an amenable blow-up rate of order 1 on Qt(S) if there exists an amenable 
blow-up rate of order 1, D(x, t), on Qt(S), so that, for some C* > 0, A > 1, and a scalar function a : [0, T] -4 

D(x, t) < \ f(x, t)\ < C*D(x, t) a ^ for all (x, t) e Qt(S). 

The condition (|6]a.i) will prove crucial in establishing our estimates in Section [3] At face value, however, it is not 
clearly motivated. To address this we specify an expansive class of functions which are simultaneously amenable 
blow-up rates of order and are reasonable blow-up scenarios given what is known about the structure of possible 
singularities in weak solutions of 3D NSE. 

Definition 7. Use the notations of Definition [6] and fix d € N. The function D(x, t) : flr{S) — > [0, 00) is an algebraic 
blow-up rate of degree d on Qt(S) if there exist functions r : (0, T] -> [0, 00), a : (0, T] -> (0, 00), and p(x, ■) : (0, T] -> 
M.[x\, X2, X3) so that, 

^"=( |p(M)l + T (J "'°" qT(5) ' 

and we additionally have, 

i. for all t E (0, T], p(x, t) is a polynomial of degree less than or equal to d and the zeros of p(x, T) are contained 
in S, 

ii. there exists A > 1 so that a takes values in [A -1 , A], 

iii. r(T) = and r(t) > for t £ (0, T), 

iv. the condition (|6]a.ii.) is satisfied. 
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Remark 8. We have presented these algebraic blow-up rates because they constitute a concrete class of amenable 
blow-up rates of order 0. To verify this, we check that condition |(6]a.i.) is satisfied, which is clear if we first expand 

log D(x,t) as, 

log(D(x, t)) = -a(t) log (\p{x, t)\ + r(t)) , 

and, additionally observe that, 

max{log \p(x, t)\, log r(t)} < log(|p(x, t)\ + r(t)) < log(2) + max{log \p(x, t)\, log r(i)}. 

Recalling the fact that, if f,g G SMO, then || max{/,g}|| BM o < 2(||/|| B mo + HsIIbmo) (BMO is a lattice), we are 
able to conclude by applying the implication (TT) in conjunction with Lemma [5] 

Our main lemma contextualizes amenable blow-up rates to their application in Section [3] As presently defined, the 
blow-up rates require that f(x,t) (as given in Definition [6} is bounded away from zero. To accommodate functions 
possibly not bounded away from zero we employ an auxiliary function, F = \J e m + f 2 , where e m is determined by 
the order of the amenable blow-up rate in question. Based on the sign of / we define, 

F _ [F, if / is non-negative, p _ fe m . if / is non-negative, 
[e m , otherwise, \F, otherwise. 

The factorization F sl£n f = F + (F-)~ 1 will allow us to independently impose blow-up assumptions on a function's 
positive and negative components. 

Lemma 9. Fix R > 0. Let S' c K 3 be a set of Lebesque measure zero for which S' n B(0,R o ) C B(0,R' ) where 
< R' Q < R Q . Suppose the function f : fi T (S") -> R satisfies f € L°°((0, T};L 1 (B(0, R a ))). Fix (f> G C£{B(0, R )) and let 
Sbea set of measure zero so that S' n B(0, R ) C S c 5(0, i? )- 

a. If J\q t (s) exhibits an amenable blow-up of order m Qt(S), then, 

sup ||^(x)log|/(x,t)||L MO < oo. 

0<t<T 

b. /fF + |g T (5) and F_\q t ( S -) exhibit amenable blow-up rates of order in Qt{S), then, 



sup 

0<t<T 



< oo. 

BMO 



c. Jjf F + |q t ( S ) and F_\q t ( S } exhibit amenable blow-up rates of order 1 in Qt{S), then, 



sup 

0<t<T 



^ll^loglogFOM) 



< oo. 

SA/O 



Proof. (Part a.) Let D(x,t) be the amenable blow-up rate exhibited by /|q t (s). We will freely reference the constants 
associated with D in the statement of Definition [6] and subsequently make the abbreviations Qt{S) = Qt(0, Ro, S) 
and B = B(Q, R ). Define an extension, F, of /|q t(S ) to fl T {S) by F = |/| on Q T (S) and F = 1 on B c x (0, T]. Set 
D(x,t) = max{D(x, t), M^ 1 } and C» = max{C*, M 2 }. These definitions ensure that, for all G Ot(5')/ 

-i-5(a;,<) < F(x,t) < C*D(x,t), 

and, taking logarithms, 

log (D(x, t)) - log < log (F(x, t)) < log (D(x, t)) + log C* . 
Recalling the implication ||TJ, we see that, 

\\logF{x,t)\\ BM o < C(M 0) C,). 
Because | / 1 = P 1 on the support of <fi we have, 

H0 lo g|/l 1 1 BMO = H^log-PHsAfO 

1 



<C||V0|U ||logP||sA/o + 



log F dx 



\B\ 

where we have applied Lemma |4] The average appearing above is uniformly bounded in time; this follows from 
the facts that A/f 1 < |/| and / G L°°((0, T]; P 1 (i?(0, Ro))). 
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(Part b.) We begin by remarking that, for all i/gl, 

(signy) log(Vl + ?/ 2 ) - 1 < , ^ „ log(\/l + y 2 ) < (signy) log \A + 2/ 2 ) + 1. 

V 1 + J/ 2 

Taking y = |/(a;, f ) |, in light of implication ((TJ it is sufficient to work with the function (sign /) log F, where F, and 
also F + and F-, are defined as in the comments preceding the statement of the lemma (taking m = 0). Those 
definitions were motivated by the fact that F slgn f = F + (F- ) ~ 1 and, so, 

log (F si ^f) = log(F + )-log(F_). 

Let G be the extension by 1 of (fF^ 1 log(F))|Q T (s) from Qt(S) to ^It(S) and adopt the notation from the proof of 
part a. of this lemma for definitions of functions analogous to those mentioned above where we designate by D± 
the blow-up rates exhibited by F±. Then, on Qt(S), we have, 



and, therefore, 



log(F+) - log(F_) - 1 < G < log(F+) - log(F_) + 1, 



\\G\\bmo < ||log(F + )|| BM o + ||log(F_)|| B MO+2. 



Appealing to part a. of the lemma, the dominating quantity is itself bounded uniformly in time. Extending this to 
the desired estimate proceeds directly, 



tjlog(F) 



<C||V0|U ||G|| 



BMO 



\BMO 



— J log Gdx 



which, again, is finite by assumptions on /. 



(Part c.) We begin similarly, noting, 



(sign /) log log F - 1 < — log log F < (sign /) log log F + 1. 



This leads us to consider, 



log (logF 



sign/ 



= log (logF), + lo § ( l °S F 



where, 



(logF) =max<{ (logF) sign/ ,l }> and (logF)_ = min I (logF) SIgnj , 1 



\ si S n / 



Recalling that F± are defined with m = 1, it is easy to check that (log F)± = \og(F±), and, therefore, 

/ 



F 



log log F = loglog(F + ) — loglog(F_ 



A further observation is that, taking D±(x, t) to be the amenable blow-up rates exhibited by F±, and setting k± = 
log(A±Mx,± + log C»,±) (these are constants associated with D± as in Definition [6] a.), we have for (x, t) G Qt, 

loglog(C^±(x,i) Q±(t) ) < \og\og{D ± {x,t)) + k* ± , 
and, recalling the comparability condition from Definition [6]c, this implies that, 

logIog(D±(x,t)) <F ± < logIog(D±(a:,*)) + k* ± . 
The conclusion now follows in the same fashion laid out in part b.; we omit the details. □ 

3. Main Results 



Weak solutions to the 3D incompressible Navier-Stokes equations are functions which satisfy (distributionally) the 
following system of PDEs, 



(3) 



dtu + (u ■ V)u = — Vp + vAu, 
V ■ u = 0; u(x, 0) = u e H, 
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where H is the L 2 closure of the divergence free test functions and the initial datum is understood in the sense of 
weak continuity (see J7) for a precise definition). The evolution of the vorticity, oj = V x u, is of special interest to 
us and satisfies, 



(4) 



d t uj + (u ■ V)w = (w ■ V) + ^ Aw, 
V-w = 0. 



For simplicity, consider a weak solutions on R 3 x [0, oo), emanating from the initial data uo that - in addition - 
satisfy the following: uj = V x u a G L 1 n £ 2 . 

The standard regularity results for weak solutions (cf. |7]) are consequences of the energy inequality and consist of 
the a priori bounds, 

SU P IMllw) < oo and / ||Vu|| i2(M 3) dt < oo, 

0<t<T> Jo 

for any T' > 0. In addition, since w G L 1 , we have (J6)), 

sup ||w|| L i (R 3) < oo. 

0<t<T' 

Note that since u; G £ 2 , our weak solution locally-in-time coincides with the smooth solution; let T be the first 
(possible) singular time. 

Choose (and fix) R > - a 'macro-scale' - with the property that B(0, R ) x {T} intersection the singular set is 
nonempty. Fix < e < i?o- Our estimates are intended for integrals over the spatial set B(0, i?o — e) and localization 
is achieved via multiplying equation ||4) by a smooth cut-off function ip satisfying, 

IWI c 

suppi/> C B(0,Rq), ^ = lonB(0,Ro-e), - — — < - for some p G (0, 1) and < ib < 1. 

V' p e 

Instead of studying the evolution of \uj\ and \uk\ directly we introduce an auxiliary function. This approach is an 
adaptation of that taken by Constantin in [6]. Define q(y) = yl + y 2 : M. —> R and let Wk = q(u>k) '■ ^ From 
these definitions it is immediate that, 

(5) |wjb| < w k , -1 < q'(cJk) = — < 1, and < g"(w fc ) = — < 1. 

Before proceeding to the main results we recall several elementary properties of the logarithm and its deriva- 
tives, 

(6) 0<-log"(j/)<log'(y)<l, ylog'2/<l, and < log'(y) + ylog"(y). 

To allow uik to have varying sign we introduce the functions Wk,+ and W).,- in the same manner as 10. For conve- 
nience we recall the notation A t (y) = {x : \u)(x, t)\ > y}. 

We include two theorems - one each for the case of solutions possessing vorticity components which exhibit 
amenable blow-up rates of order and of order 1 - and their proofs. The proof of the order case adapts directly to 
the case of blow-ups of order 1 and the proof of the second result is accordingly terse. 

Theorem 10. Let ubea Leray solution to @ oh R 3 x [0, oo), and suppose that - in addition -wo = Vxu eL 1 n L 2 ; let 
T be the first singular time and S the singular set ofuj at time T. Fix positive values Rq and e so that ^ S PI B(0,Rq) c 
B(0, R —e)-i.e. there are singular points in B(0, R ) but all such points are in B(0, R —e)-and let Qt(S) — Qt(0, Ro,S). 

(i) If Wk,+ \Q T (s) an d Wk,-\Q T (S) ea ch exhibit amenable blow-up rates of order in Qt(S), then there exists a positive 
value Mfe so that, 



sup / \u h (x,t)\ log ( + \u k {x,t)\ 2 ) dx < M k . 

te[0,T] JB(0,R -e) 

(ii) If for all k, the premise of part (i) is satisfied, then there exists a positive value M so that, for t G [0, T], 
1 V\ . M 



y ° Z ( At (^ l|wWlli ~( s (° A - e »)) 



Ki)||L°o( B (0,.Ro-e)) log (l + \\u{t)\\ L ~, {B{0 ,Ro~e))) 



8 Z. BRADSHAW AND Z. GRUJIC 

Remark 11. As will be seen, our energy inequality-type method depends heavily on classical techniques and we 
thus need to consider smooth solutions. It is possible to extend the estimates to some weak solutions by considering 
a sequence of smooth approximations - a good choice to use here is the method of retarded mollifiers given for 
suitable weak solutions in 111 as it allows one to recover from a velocity-level approximation scheme information 
about the vorticity see [6| for more discussion - but, as our estimates hinge on assumptions beyond the initial data, 
Mo, we would have to ensure these are met by a convergent sub-sequence of approximate solutions. 

Proof. In virtue of our smoothness assumption, the evolution of ipw k log(wfe) can be established from the evolution 
of Wfc by first writing, 

d t (ipw k log w k ) 

= i'q'(ui k )( log w k + w k log' w k ) (d t uj k ) 

= ipq' {w k )(k>gWk + w k log' w k ){y/\u k - (u ■ V)w fe + (uj ■ W)u k ), 
and then deriving (noting tacit summation over terms involving indices other than k), 

d t (i/jw k log w k ) - vipAw k ( log w k + w k log' w k ) + H'q"(^k)(di^ k ) 2 ( log w k + w k log' w k ) 

= il>(uo ■ V)ufc</(<jfc)(log Wk + Wk log' Wk) - ip(u ■ \7)wk(logw k + Wk log' Wk). 

By integrating in space and time and dropping the positive quantity involving q" from the left hand side we obtain 
the following energy inequality-type estimate, 

(7) J tp(t)wh(t)logWk(t) dx - v J J ipAw k (\og Wk + w k log' w k ) dx ds 

J J ij)(oj ■ \J)u k q'(uJk)(\ogWk +w k \og' w k ) dxds 
i/j(u ■ V)w k ( log Wk + Wk log' Wk) dx ds 



< 



+ J ^ Wkfi log w k fi dx. 

The properties in |(6) enable several key cancelations. For the dissipative terms, integration by parts reveals that, 

J j \d?w k )ip log w k dx ds = J J \w k (dfa) log w k + w k (djip)(djW k ) log' w k ^j dxds 

-J J {djWkf iplog' w k dx ds, 

and, 



(djW k )tpw k log' w k dx ds = - j j w k {djil))(djWk) log' w k dx ds 

J (djWk) 2 -^ ( log' w k + Wk log" Wk^j dx ds. 



A cancelation occurs upon adding the above equations leaving us with, 

(d 2 w k )ip I log Wk + w k log' Wk ) dx ds 



= v J J (djWk) 2 ^ \2log' w k + w k \og" Wk^J dx ds - v J J w k (d 2 ip) log w k dx ds. 

Again noting |(6), the integrand of the first term is positive and, therefore, can be dropped from the left hand side 
of estimate 0. The second term can be dominated by an a priori finite quantity arising from the standard energy 
inequality for weak solutions in conjunction with the fact that our spatial integral is over a set of finite measure. 
More precisely, 

(8) / / w k (d ftp) log w k dx ds < C / w\ dx ds 

Jo J Jo Jb(o,r ) 



cf [ \uj k \ 2 dxds + T\B(Q,R ) 
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The integrals arising from the transport term in © also enjoy a substantial cancellation, 

J j ip u ■ Vwk ( log Wk + Wk log' Wk*J dx ds 

ip Uj (djWk)wk log' Wk dx ds — J J ip UjWk (djWk) log' Wk dx ds 
UjWk log uikdjip dx ds 
UjWk log Wkdjip dx ds. 



Noting that log(y) < Ay 1 / 4 , 



< C 



UjWk logWk\djijj\ dx ds 



<C \\V /4 u\U\\w k \\Li 4 (B(0 Rn ^WwkWm 



<C sup ||«;*||i( 4 (B(0)Jlo)) 

0<t<T 



L 1 (B(0,-Ro)) " u/ " I l^ z (B(0,-Ro)) 
T 

||V(^ p/4 M)|| 2 ||w fc || i 2 (B(0ii?o)) , 



where we have used Holder's inequality and the Sobolev inequality. From here, a commutator estimate on the 
gradient of the localized velocity as well as reasoning along the lines given for ((HJl allow the extension of the above 
estimate to one in terms of a priori finite quantities. 

The integral on the last line of the right hand side of is finite by our assumptions on the initial data. 
At this point, based on 10, we have established, 

J(tpWk log Wk){t) dx < J J tp u> ■ Vufc q'(w k ) ^logw fc + w k log' w k ^j dx ds + R, 

where R is comprised of those a priori bounded quantities accumulated in the preceding estimates. Noting that, 

i>q'(uik)w k log' w k < 1, 

an a priori bound follows for that part of the as of yet unbounded quantity leaving us with, 



(i/jw k log w k )(t) dx < 



UJk 

ip uj ■ Vitfe — log Wk dx ds 



R. 



The T-0--BMO duality and the Div-Curl lemma justify the following chain of inequalities, 



uj ■ Vuk — log w k dx ds 



< / \\ui ■ Vitfe 



H 1 



UJk 

-tjj — log w k 

Wk 



< [ sup 

,0<t<T 



< ( sup 

,0<t<T 



UJk 

i> — logwfc 

Wk 



Ip logltffc 

Wk 



BMO 



ds 

BMO 
T 



/ ||w-Vufc|| H1 ds 
Jo 



T \ 1/2 / r T 

2 





ML dt 



Vu\\l dt 



1/2 



BMO/ / 

By Lemma|9]and the standard regularity of Leray weak solutions all of the above are finite and we have thus estab- 
lished that, for all < t < T, \\ipWk log Wfc)(£)||Li(R3) is majorized by time-independent a priori bounded quantities, 
the sum of which we label Mk, completing our proof of part (i) of the theorem. 

The second conclusion - part (ii) of the theorem - follows in two steps. The first is just Chebychev's inequality, 



Vol(A t f i||w(t)|| L co (fl(0)Ho _ e)) ] J < 



cJo{t) 



|w(i)||i~( S (0,i{o-6)) log (1 + ||w(i)||i=c (1J( o >iJo _e))) ' 
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where, 
for, 
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I (t) = I Hlog(l + |w(a;,t)|) dx, 

I {\u(x,t)\>\(t)}nB(0,Ro-e) 



Ht) = ■^\\u(t)\\L°°(B(0,R o -c))- 



The second step is largely technical and reconciles the difference between the vorticity-component determined 
quantities above and the vectorial magnitudes apparent in Io(t) in order to bound Io(t) in terms of the quantities 
Mfc. Tacitly summing over i, j we have, 

Io(t) <C \uj\iogJu^ dx. 

JB(0,R o -e) 

An explicit reduction illustrates the argument (for simplicity we take j = 1 and integrals over the indicated sets 
intersected with 5(0, Rq — e)), 



\u>i \ log + cu 2 + lo 2 dx 

(|w 2 | + I^Dlog \j%u^2 dx + / wi log \J 2w\ dx. 

Applying the same reasoning to the first integral above and then repeating for all values of j eventually yields, 

Io(t) < C J Wi log ^Jcwi dx < C J Wi dx + C J Wi log Wi dx < CMi. 



□ 



The energy inequality-type construction used to prove the previous theorem also works if we substitute log" 1 u)k 
(as defined in Section [2} in place of log Wk ■ This allows the application of amenable blow-up rates of order 1 in 
conjunction with Lemma[9]c. To ensure things are meaningful we modify our definition of q so that q{y) = \fe + y 2 
(so, now, Wk = y/e + uify and refer to l(2j to define Wk,+ and Wk.- taking rn = 1. 

Theorem 12. Let ubea Leray solution to ||3) on M 3 x [0, oo), and suppose that - in addition -u = Vxu e L 1 n L 2 ; let T be 
the first singular time and S the singular set at time T. Fix positive values Rq and e so that ^ S (~l B(0, R ) c B(0, R — e) 
- i.e. there are singular points in B(0, R ) but all such points are in B(0, R — e) - and let Qt(S) = Qt(0, Ro, S). 

(i) If Wk,+ \Q T (s) and Wk,-\Q T (s) ea ch exhibit amenable blow-up rates of order 1 in Qt{S), then there exists a positive 
value Mk so that, 



sup / \uk(x, t)\ loglog (^e+ \uj k (x,t)\ 2 ) dx < M k . 

t£[0,T] JB(0,R -e) 

(ii) If, for all k, the premise of part (i) is satisfied, then, there exists a positive values M so that, for t e [0, T], 

M 



voiUt(^\\w{t)\\ L ^ B(0iBo - e ))y\ < 



l^(*)l|i-(S(0,Ro-e)) loglog (e+ \\u>(t)\\ L ^ [B{0<Ro _ e)) ) 



Proof. Multiplying equation ||4} by tpq'(wk)( loglogWfc + Wfe(loglog) / 'u;fe) we obtain the evolution of ipwu log log w k - 
The point-wise estimates in (|6} adapt directly to the function loglog(y) (indeed, they adapt to any number of such 
self-compositions of the logarithm) and, after integrating in space and time, all of the estimates and cancelations 
from the previous proof - except those involving the vortex stretching term - can be duplicated directly. We thus 
obtain, denoting by R some a priori finite quantity, 



(9) 



(ipWk \og\ogWk)(t) dx < 



ip u ■ Vitfc — log log Wk dx ds 

w k 



R, 



and, using Lemma|9]c., we are able to pull out of the integral the uniformly-in-time bounded quantity, 



sup 
t 



OJk 

ip — log log w k 

Wk 



BMO 



and conclude exactly as in the previous proof. 



□ 
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